Spectral properties of Bunimovich mushroom billiards 
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Properties of a quantum mushroom billiard in the form of a superconducting microwave resonator 
have been investigated. They reveal unexpected nonuniversal features such as, e.g., a supershell 
effect in the level density and a dip in the nearest-neighbor spacing distribution. Theoretical pre- 
dictions for the quantum properties of mixed systems rely on the sharp separability of phase space 
- an unusual property met by mushroom billiards. We however find deviations which are ascribed 
to the presence of dynamic tunneling. 
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Billiards play a central role in the investigation of sys- 
tems with regular, chaotic, or mixed dynamics When 
quantum chaos established itself as a new field, billiards 
were used as a paradigm for both theoretical and ex- 
perimental research [2[- The billiard considered in this 
Rapid Communication is from the family of mushroom 
billiards suggested by Bunimovich [3[ as a generalization 
of the stadium billiard Compared to conventional 
mixed systems they have the particular property that 
their phase space is sharply divided into one regular is- 
land and the chaotic sea, whereas usually the islands of 
regularity are typically surrounded by a layer of infinitely 
many islands. In the simplest case they consist of a semi- 
circular hat with a symmetrically attached rectangular 
stem as shown in Fig. QJa). All regular orbits of mush- 
room billiards are orbits of the semicircle billiard with a 
conserved angular momentum which stay in the hat for- 
ever. Orbits of particles with the same angular momen- 
tum form a semicircular caustic. There is one critical 
caustic, which rigorously separates the orbits into regu- 
lar and chaotic ones; its radius r c equals half the width 
of the stem [dotted line in Fig. [Ha)]. Particles moving 
in the hat with a larger or equal caustic stay there for- 
ever, whereas those with a smaller one eventually enter 
the stem and therefore are chaotic. To our knowledge 
this clear separation of the phase space has been found 
before only for classical maps [5]. Due to this unusual 
feature of the phase space classical mushroom billiards 
are of interest with respect to different aspects [6| . How 
comes about the separability of the classical phase space 
manifest in properties of the corresponding quantum bil- 
liard? This question motivated us to investigate a quan- 
tum mushroom billiard experimentally using the analogy 
between quantum and microwave billiards 0, 0] . 

The geometry of the quantum billiard we investigated 
is shown in Fig. [2(b). To avoid the effects induced by 
the superposition of two parity classes, we used a desym- 
metrized mushroom billiard - i.e., one with a quarter- 
circle for the hat. Moreover, the stem is chosen triangular 
instead of rectangular in order to eliminate bouncing ball 
orbits as seen, e.g., in the stadium billiard [8]. The stem 
width is two thirds of the radius R of the hat, and its in- 



FIG. 1: (Color online) (a) Typical shape of a mushroom bil- 
liard with the critical caustic shown as a dotted line, (b) 
desymmetrized version (quarter-circle as hat) with triangular 
stem, and (c) photograph of the corresponding experimental 
microwave resonator. The positions of the six antennae in the 
experiment are indicated in (b). 



ner angle equals 45°. We verified with a rigorous analysis 
[3[ that the classical phase space is still sharply divided 
into a regular and a chaotic part. The invariant measure 
q c of the chaotic part is 82.9 % of the phase space vol- 
ume. The eigenvalues of the quantum mushroom billiard 
were measured with a flat, cylindric microwave cavity of 
lead plated copper [Fig. [2(c)] of 5 mm height. The hat 
has a radius of R = 0.24 m. In order to obtain a large 
and reliable set of resonance frequencies we performed 
the measurements at 4.2 K, where the cavity is super- 
conducting [1, Using a vectorial network analyzer 
(VNA) we collected complex transmission spectra with 
six different antennae up to 22 GHz with a sample rate 
of 100 kHz. The antenna positions are distributed over 
the whole area of the billiard [Fig. [2(b)]. Figure [2] shows a 
part of the measured transmission spectra for five differ- 
ent antenna combinations. The spectra have the unusual 
property that they exhibit sequences of resonances, which 
are separated by large gaps with no resonances (arrows in 
Fig. [2]). This bunching effect indicates that there are two 
well distinguishable frequency scales. It will be discussed 
in detail below. 

A set of 938 resonance frequencies was obtained from 
the spectra in agreement with the expectation from 
Weyl's formula [10[. On the basis of this large data set 
we first considered conventional statistics and compared 
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FIG. 2: Part of five transmission spectra between 9 and 
10 GHz, measured with the microwave resonator shown in 
Fig. [2(c). The curves labeled with \Sij\, j = 2,..., 6 are ob- 
tained with the antennae 1 and j whose positions are indi- 
cated in Fig. \Hb) . The spectra have been displaced from each 
other along the y axis for illustration and are plotted in a log- 
arithmic scale. One notices a bunching behavior - i.e., large 
gaps marked by arrows - between stretches of resonances. 



them to theories applicable to the quantum spectra of 
mixed systems, such as, e.g., the Berry- Robnik statistics 
111 ]. First we computed the number of levels N(f) - i.e. 
the integrated resonance density p(f) - and determined 
its smooth part iVweyi(/)- The resulting fluctuating part 
AW(/) = N(f) - N Weyi (f) is shown in Fig. The 
most striking feature is a beating pattern caused by the 
superposition of two oscillations, which in fact reflects 
the bunching effect observed in Fig. [2j In the stadium 
billiard Q oscillations are caused by the so-called bounc- 
ing ball orbits. To find the origin of the beating observed 
in Fig. 03 we computed the length spectrum |pfl uc (x)| 
shown in Fig. UK a), which has peaks at lengths x of pe- 
riodic orbits. Here, pnuc(x) is the Fourier transform of 
the fluctuating part of the resonance density pnuc(k) as 
a function of wave number k = 2tt/ / c, where c is the ve- 
locity of light. In order to get information on the nature 
of the periodic orbits related to the peaks in the length 
spectrum, we also computed the latter for the regular 
orbits in the hat [8]. For this purpose we first deter- 
mined the eigenvalues and eigenfunctions of a quarter- 
circle billiard of the same radius as the hat. These are 
indexed by a radial quantum number n and an angular 
momentum one m. The length spectrum in Fig. 2Kb) 
has been obtained by considering only eigenvalues with 
eigenfunctions, which are localized between the critical 
caustic [Fig. [2(a)] and the circular boundary - that is, 
eigenfunctions with a sufficiently large m. By compari- 
son of Fig. [41(a) with Fig. 2Kb) we see that all peaks of 
regular orbits in the hat of the mushroom billiard are re- 
produced. The remaining peaks in the length spectrum 
are thus associated with chaotic orbits. Around 0.7 m 
we observe a pair of closely lying dominant peaks. They 
correspond to the shortest regular periodic orbits. Af- 
ter subtracting their contribution to iVfl uc (/) the beating 
vanished. This cause of beatings in circular or spherical 
geometries by short dominant orbits of approximately the 
same length has already been observed in metal clusters 
[12[ and termed the supershell effect. It has also been 
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FIG. 3: Fluctuating part of the number of resonances below 
a given frequency /. The curve oscillates and shows a beating 
behavior with some noise at the nodes marked by arrows. 
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FIG. 4: (a) Experimental length spectrum of periodic orbits 
of the mushroom billiard. The two peaks of lengths of about 
0.7 m are due to two regular orbits in the hat of the mush- 
room, and that at length 1.12 m is due to a chaotic one. (b) 
The computed length spectrum of a quarter-circle billiard is 
shown flipped upside down, where only eigenvalues were con- 
sidered which correspond to angular momenta larger than the 
critical one. 



stressed in nuclei [13], nanowires [141 ]. and recently in 
trapped dilute Fermi gases [l5[ . In mushroom billiards it 
is induced by regular orbits coexisting with chaotic ones 
whose measure is determined by the depth of the stem. 

Next we consider the nearest-neighbor spacing distri- 
bution (NND). Figure 03(a) shows the experimental NND 
with that for Poisson statistics and a Gaussian orthogo- 
nal ensemble (GOE) describing regular and chaotic sys- 
tems, respectively. As to be expected the experimental 
NND coincides with neither of them. A comparison of 
the NND with the one derived by Berry and Robnik 
for mixed systems yields good agreement for q c = 82.9 % 
as given above. Yet, even though the requirements for 
its applicability are best met with mushroom billiards, 
most surprisingly the agreement is not better than for 
other mixed systems. Especially, for small spacings s de- 
viations are visible in Fig. 03(a) . This can be attributed 
to dynamic tunneling between classically se par ated re- 
gions of phase space in the quantum regime [lg, LL7|, LL8| . 
In addition, we observe a peculiarity in the shape of 
the NND. Namely, the experimental histogram result- 
ing from a fairly small binwidth exhibits a profound dip 
at s = 0.7. This dip is statistically significant, and it 
is present also in NNDs derived from different subsec- 
tions of the spectrum. When including the contribution 
of the pair of dominant short periodic orbits to N^ uc (f) 
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FIG. 5: (a) Nearest- neighbor spacing distribution of the 
mushroom billiard together with those for Poisson statistics 
and the GOE (solid lines) expected for systems with regular 
and chaotic dynamics, respectively. Note the dip at s = 0.7. 
It vanishes (b) when the contribution of the dominant peri- 
odic orbits in Fig. [4] to the resonance density is subtracted. 
The dashed line in both (a) and (b) shows the Berry-Robnik 
prediction for a mixed system with q c = 82.9 % chaos. 



in the usual spectral unfolding procedure [8], the dip is 
no longer present [Fig. 03(b)]. It should be noted that 
this result is contrary to the general belief that only long 
periodic orbits control short range spectral properties of 
a quantum system. 

While in the semiclassical limit the eigenvalue spec- 
trum of a mixed system consists of regular eigenvalues 
with eigenfunctions localized on the regular islands - 
i.e., for the desymmetrized mushroom billiard those of a 
quarter-circle billiard with eigenfunctions localized in the 
hat - and of chaotic ones with eigenfunctions distributed 
over the whole billiard [l9[, this must not be true any- 
16 , [ItJ in the quantum limit. However, since for 
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mushroom billiards there is no layer of islands between 
the regular and chaotic parts of the phase space, their 
eigenstates should be classifiable as regular or chaotic 
[la l20| even in the quantum regime. In order to investi- 
gate this further, we measured the electric field strength 
intensity - i.e., squared eigenfunctions of an enlarged 
copy of the mushroom billiard - for several resonance fre- 
quencies using the perturbing bead method at room tem- 
perature [21]. An example for a squared regular eigen- 



function is shown in Fig. [61(a). It is very similar to that 
of the quarter-circle with n = 3 and m = 44. A chaotic 
eigenfunction is plotted in Fig. [51(b). However, for such 
eigenfunctions we even find traces of regularity in the 
field pattern in the hat, and for the regular eigenfunc- 
tions the intensity in the stem is nonvanishing. This indi- 
cates that there is a dynamic tunneling [18[ , such that the 
classification into regular and chaotic eigenstates is only 
asymptotically correct. Indeed, we also observed some 
rare mixed eigenfunctions whose intensity is equally dis- 
tributed over the whole billiard area, while the field pat- 
tern in the hat is very similar to that of a regular eigen- 
function [Fig. EKc)]. In Fig. [51(d) the averaged intensity 
distribution of 239 chaotic eigenfunctions is shown. Its 
classical counterpart is the probability Pq to find chaotic 
orbits at a certain position in the billiard. It is constant 
as in completely chaotic systems in that part of the bil- 
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FIG. 6: (Color online) Measured intensity distributions. The 
eigenfunctions are usually either regular (a) or chaotic (b) and 
correspond to eigenfunctions of a quarter-circle billiard in the 
former case. We also found some mixed eigenfunctions like 
the one in (c). In (d) the averaged distribution of 239 chaotic 
eigenfunctions is shown. The radial intensity profile in the 
hat is plotted in (e) between the quarter- circle center and 
the circular boundary (points), and compared to the classical 
probability of finding chaotic orbits (solid line) . The decrease 
beyond the critical caustic at r/R — 2/3 is revealed in the 
measurements, and deviations at the boundaries and at the 
knee are due to quantum effects. 



Hard which is accessible to chaotic orbits only - i.e., in 
the stem and in the hat for < r < r c , where r is the 
distance from the circle center and r c is the radius of the 
critical caustic. Interestingly, for r c < r < R it decreases 
as arcsin(r c /r) [3], which is a consequence of equipartion 
in the mixed phase space. The radial profile of Pc is com- 
pared to the experimental results in (e) - revealing the 
decrease of intensity in the hat. The abrupt change of the 
distribution when crossing the critical caustic is smeared 
out due to the finite wavelength. At the boundary of the 
billiard quantum effects [22[ also lead to deviations from 
the classical behavior. 

Finally, for a fixed n, the eigenvalues of the quarter- 
circle become asymptotically equidistant for large fre- 
quencies. Thus regular modes in the superconducting 
mushroom billiard were obtained by identifying chains 
of equidistant eigenvalues in the measured spectrum and 
comparing them to the computed ones of the correspond- 
ing circle billiard. We found regular (chaotic) periodic or- 
bits in the length spectrum of the regular (chaotic) eigen- 
values and, though strongly suppressed, also chaotic (reg- 
ular) ones appear. This shows again that there is an in- 
teraction via dynamic tunneling between the regular and 
chaotic parts of the classical phase space. The deviations 
of the measured regular eigenvalues of the mushroom bil- 
liard from the computed eigenvalues of the quarter-circle 
billiard shown in Fig. [71 quantify the dynamic tunneling. 
Most interestingly they vary unidirectional. Indeed, for 
all families of equal radial quantum number n the devi- 
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ing process. 
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FIG. 7: Deviations of the experimental eigenfrequencies of 
the mushroom billiard from those of a quarter-circle billiard 
versus frequency. The symbols (points, triangles, squares) 
refer to regular states. Eigenvalues with equal radial quan- 
tum number n are connected by lines. For a given n there are 
almost no deviations for high frequencies where the eigenfunc- 
tions are localized close to the circular boundary. For smaller 
frequencies, the influence of the stem leads to increasing de- 
viations pointing to the importance of dynamic tunneling. 

ations are close to zero for large frequencies - i.e., large 
values of m - where the eigenfunctions of the quarter- 
circle billiard are localized close to the circular bound- 
ary and therefore almost not influenced by the stem. 
When approaching the frequency corresponding to the 
critical angular momentum from above, the deviations 
increase as then the eigenfunctions are localized closer 
and closer to the critical caustic and the field distribu- 
tion is increasingly distorted towards the stem as can be 
seen in the intensity plot in Fig. [6ja). This is exactly 
the dynamic tunneling. There is no general quantitative 
explanation for the measured deviations yet. Our results 
should, however, should inspire forthcoming theoretical 
considerations of dynamic tunneling. 

In summary we have shown that the spectral proper- 
ties of mushroom billiards are significantly affected by the 
shortest regular orbits in the hat. They cause a super- 
shell structure in the level density. As compared to other 
systems exhibiting the supershell effect, in mushroom bil- 
liards the degree of chaoticity can be tuned by varying 
the depth of the stem, such that they allow the study of 
supershell effects for an arbitrary degree of chaos. Sur- 
prisingly, the shortest (not long) periodic orbits lead to 
a substructure in the NND. The eigenstates of the mush- 
room billiard may be separated into regular, chaotic, or 
- though rare - mixed ones. That the latter are rare and 
the behavior of the averaged intensity distribution of the 
chaotic eigenfunctions in the mushroom hat are manifes- 
tations of the separability of the classical phase space in 
the spectral properties of the quantum billiard. Still, dy- 
namic tunneling is present and can be observed, e.g., in 
the field distributions, in the spectral properties, and in 
the deviations of the regular eigenvalues from those of the 
corresponding circle billiard. Since the tunneling barrier 
is of a simple structure for mushroom billiards, they are 
convenient systems for the study of the dynamic tunnel- 
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